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Abstract. Gravitational radiation from the galactic population of white dwarf 
binaries is expected to produce a background signal in the LISA frequency band. 
At frequencies below 1 mHz, this signal is expected to be confusion-limited and 
has been approximated as gaussian noise. At frequencies above about 5 mHz, the 
signal will consist of separable individual sources. We have produced a simulation 
of the LISA data stream from a population of 90k galactic binaries in the frequency 
range between 1-5 mHz. This signal is compared with the simulated signal from 
globular cluster populations of binaries. Notable features of the simulation as well 
as potential data analysis schemes for extracting information are presented. 



1. Introduction 

The Laser Interferometer Space Antenna (LISA) is a planned space-based gravitational 
wave detector. It will be sensitive to gravitational radiation in the frequency band 
from 10^^ — 1 Hz, with sources ranging from supermassive black hole binaries to 
stellar mass compact binaries. One of the most abundant sources will be the galactic 
population of close white dwarf binaries. These sources will dominate the LISA signal 
from 0.1 mHz up to about 3 mHz and are expected to produce a confusion limited 
background below 1 mHz. The signal will transition to individually resolved sources 
above ~ 5 mHz. Data analysis strategies to accommodate the close white dwarf 
binaries vary according to the nature of the signal. Below 1 niHz, the signal will 
be isotropic and effectively gaussian with > 100 sources per frequency bin for a one 
year observation. Above 5 mHz, the signal will individual sources widely separated 
in frequency. The transition region between 1 and 5 mHz presents unique challenges 
for data analysis. In addition, sources of gravitational radiation in globular clusters 
have the greatest angular resolution in this transition zone. It is therefore of use to 
develop a simulation of the LISA data stream from galactic white dwarf binaries in 
the transition zone between 1 and 5 mHz test data analysis techniques for handling 
the foreground signal while identifying other signals of interest such as globular cluster 
binaries or high mass ratio inspirals into supermassive black holes. We have generated 
just such a simulation with 90,000 white dwarf binaies and describe in Section [3 the 
population synthesis and in Section|31the data stream simulation. The simulated data 
stream from select globular clusters is described in Section 0] 
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Table 1. Probabilities of binary types at birth. 
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2. Population Synthesis 

The population of galactic white dwarf binaries is synthesized using a phenomenologi- 
cal approach. We assume a uniform birthrate over the past 10^° yr, and then assign a 
binary type at birth according the simulations of Nelemans et al (2001). The birth of 
the binary is considered to occur after the second white dwarf is formed. The binary 
types are combinations of four basic white dwarf types - Helium, hybrid Carbon Oxy- 
gen (with a thick Helium envelope), Carbon Oxygen, and Oxygen Neon. The binary 
types and their probabilities at birth are shown in Table [3 Once the binary types 
have been determined, individual masses are assigned to each white dwarf based upon 
its type. These are drawn from mass distributions derived from Iben and Tutukov 
(1986). The orbital period, Q (in hours) at birth is chosen by assuming a normal 
distribution with unit variance in log (Q) with the mean period Q related to the mass 
of the secondary using m = 0.0625 log (Q) + 0.2625 if the secondary is a Helium white 
dwarf and log (Q) = 0.4 if it is any other type. Finally, the binary is assigned a birth 
date between the present and 10^° years ago. The binary is allowed to evolve under 
the influence of gravitational radiation to determine the present orbital period (P) 
using 

P = [q'" ^\k,t^"' (1) 

with 

96 ,8/3 G^/^M,M2 
ko = — (27r) -— (2) 

5 c5(Mi + M2)'/^ 

where Mi and M2 are the masses of the white dwarfs (Peters and Mathews, 1963). 
Finally, the binary is retained as part of the synthesized population if its orbital period 
is less than 2000 s and greater than the period at Roche lobe contact. 

Once the masses and orbital period of the binary have been established, it is then 
assigned a position in the galaxy according to the probability distribution functions 
in the cylindrical coordinates R and z given by: 

P{R) = i-e-^/^" (3) 

P(z) = -^sech2(^/zo), (4) 
2zo 
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where Rq = 2500 pc and zq = 200 pc. The angular coordinate is chosen from 
a flat distribution and the coordinates are converted to heliocentric spherical polar 
coordinates, (r, 0s,<^s), with (jj^ measured from the vernal equinox and measured 
from the normal to the ecliptic. The orientation of the binary is described by the 
direction of the angular momentum vector. The angular coordinates of this vector are 
assigned using flat distribution functions and are {0l,4'l) in the same heliocentric 
coordinate system. Finally, the orbital phase of the binary ((/jq) is also chosen 
randomly. The population of binaries is then described by 9 parameters - the masses 
of the components (Mi and M2), the orbital period (P), the position in heliocentric 
coordinates {r,9s,(j)s), the orientation {9L,(j)L), and the orbital phase ((/io)- 

The number of binaries to include in the population is chosen to be 90,000 which 
is close to the number determined by using the population estimates of Hils, Bender, 
and Webbink (1990) and (following their lead) reducing this number by a factor of 10. 
The distribution of the masses of the binaries are shown in Figure and the number 
of sources per frequency bin for selected frequencies are shown in Figure |21 

3. Data Stream Simulation 

The LISA data stream simulated from the population of binaries is calculated using 
the long wavelength approximation of Cutler (1998). Cornish and Rubbo (2003) have 
shown that this low frequency approximation will suffice for frequencies below ~ 5 
mHz. We have calculated the data stream for both the hj and hu configurations 
of Cutler and including both the Doppler phase {(Pd) and the polarization phase 
{(fp). We have chosen for the initial configuration of LISA to be = and = 
which correspond to LISA starting at the vernal equinox at t = 0. The binaries are 
assumed to have been circularized and therefore emitting monochromatic gravitational 
radiation with frequency / = 2/P. The data stream is calculated for each binary 
as a double-precision strain amplitude measured every second for 32 x 10^ seconds. 
These data streams have then been combined to form a single stream for all 90,000 
binaries. The data stream for hj is shown in FigureOl It is of interest to note that the 
amplitude of the signal from the galactic binaries can be seen to vary annually as the 
antenna pattern of LISA sweeps through the plane of the galaxy, thus showing that 
the anisotropy of the binary distribution is measurable in this frequency band. 

4. Globular Cluster Simulations 

Although there have been a few simulations of the white dwarf binary population of 
globular cluster (Hurley and Shara 2003), a detailed synthesis of such a population 
is very difficult due to the additional complication of the internal dynamics of a 
globular cluster. Consequently, we have simply chosen to use binaries from our 
galactic population synthesis to simulate the signal from globular cluster binaries. 
Since the globular cluster binaries are more localized in space, the time variation of 
their signal should be different than that from galactic binaries and this is borne out in 
the simulations performed for five individual globular clusters containing 10 binaries 
each. The globular cluster populations were simulated by choosing 10 binaries from 
the population synthesis described in Section |21 with their positions chosen from a 
three-dimensional normal distribution centered on the position of the globular cluster 
with a variance given by the half-mass radius. The globular cluster data was taken 
from Harris (1996). The datastreams from these globular clusters and a combined 
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Figure 1. Distribution of the masses of the white dwarfs in solar masses. The 
upper histogram shows the masses of the first white dwarf formed while the lower 
one gives the second white dwarf formed. 



data stream are shown in Figure 01 Because the globular clusters are much more 
localized sources, their signals can be reduced to zero as it passes through a null in 
the antenna pattern of LISA. Searching for such nulls may provide a technique for 
identifying globular cluster binaries. 

5. Conclusion 

We have synthesized a population of close white dwarf binaries in the galactic disk and 
used this population to generate a simulated data stream from LISA. This simulation 
does not include the instrument noise for LISA, but noise can be included in the signal 
by simply adding the appropriately chosen noise spectrum to the signal. In this way, 
our simulated data stream can be used to compare the effect of different instrumental 
noise spectra on the ability of different data analysis techniques to handle the signal 
from the galactic population of white dwarf binaries. The simulation can also be 
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Figure 2. Number of binaries per 1 yr~^ frequency bin. The upper histogram 
shows the first 100 bins above 1 mHz while the lower one shows the first 100 bins 
above 3 mHz. 



combined with simulated signals for other sources of interest (as we have done with 
globular cluster sources), to test different data analysis schemes for extracting these 
other signals. 

We anticipate improving the simulation by including white dwarf binaries with 
orbital periods greater than 2000 s and incorporating the general response function of 
Cornish and Rubbo (2003) for binaries with frequencies higher than ~ 3 mHz. In the 
meantime, we will make the simulations freely available on CD for anyone interested 
in using it for testing data analysis schemes. 
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Figure 3. Full year data stream of hi strain amplitude. The time axis covers 
32 X lO'* s and the strain amplitude is less than 3 x 10"'^'^. 
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Figure 4. Data streams for NGC 104, NGC 5139, NGC 6397, NGC 6752, NGC 
7078, and the five globular clusters (in grey) compared with the disk signal (in 
black). The data streams from the individual globular clusters are plotted at 
difl'erent scales to show the details. The actual amplitudes range from 2 x lO"-^-^ 
to 1.1 X 10-21. 
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